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Abstract 

O 

Q\ . A recipe for the generalization of the Boltzmann equation to a quantum ki- 



netic equation is given for cases in which only level shift and broadening are 



c3 

a. 

f^j . considered, while coherence phenomena can be neglected. We also consider 

a: 

a specific problem: Hot luminescence from a quantum wire near the thresh- 



old for optical phonon emission. The problem is first discussed within the 
framework of the Boltzmann equation. After pointing out the failure of this 
description, the Boltzmann equations are generalized to a set of quantum ki- 
netic equations, which in turn are solved in order to describe the luminescence 
spectra. 



Typeset using REVTpjX 



I. INTRODUCTION 

The simplest way to describe kinetic phenomena is by the Boltzmann equation (BE). 
There are however cases when the BE cannot be used, and one has to use quantum kinetic 
equations (QKEs) that are in the simplest case equations for the one particle Green functions. 
If the Keldysh technique [0-0] is used, the one particle Green functions that participate in the 
description of kinetic processes are the retarded Green function (and its complex conjugate), 
and the "statistical" Green function, which is simply the field correlator. For free particles 
these Green functions are diagonal in the basis of the free particle states, and are described 
by the free particle energies e a , and their occupation numbers n a . These are the quantities 
that appear in the BE for the interacting particles in external fields. 

The description of kinetic phenomena by QKEs instead of BEs results in two effects. 
The first effect is state mixing - due to interactions and external fields the Green functions 
are nondiagonal in the basis of the free particle states. This effect is responsible for co- 
herence phenomena, of which a well known example is interband polarization described by 
semiconductor Bloch equations (§-§]]. 

The second effect is level shift and broadening - the Green functions are nonzero off the 
energy shell within some width defined by the interactions and the external fields. Level 
broadening effects in QKEs for phonons have been studied in Refs. (7]||. 

The effects of level broadening are of importance in two cases. The first case is that 
in which level broadening allows scattering processes that are forbidden by energy and 
momentum conservation. An example of such a case is that of electron-electron scattering 
in a 2DEG, in the presence of a quantizing magnetic field ||. If one tries to calculate the 
Auger transition rates responsible for Anti-Stokes luminescence, that was observed in Ref. 
1(J, using Fermi's golden rule, one encounters delta functions with an argument equal to 



zero, due to singularities in the density of states, see Ref. JTTJ . 

The second case in which level broadening effects should be considered is that of narrow 
energy distributions: When there are small energy scales in the problem that are narrower 



than the quantum width of the state. A simple example of such a case, that has one small 
energy scale, is that of an electron distribution, excited by a narrow band-width laser. The 
problem becomes more complex when another small energy scale is added, for example: The 
electrons are excited to the close vicinity of some threshold. A problem of this nature will 
be dealt with in this paper. 

In the time domain, level shift and broadening lead to non-Markovian effects in the 
scattering processes [112,131, 6fl. These effects are important in the case of short excitation 



pulses [PH17 



If state mixing induced by coherent excitation can be neglected, and only shift and 
broadening effects are considered, the QKEs can be presented in a form similar to the BEs. 
We will give a general recipe how to write such equations, and will use this recipe when 
considering the following problem: Hot luminescence from a quantum wire (QWR) due to 
electrons that are excited by a narrow band, noncoherent laser, to the close vicinity of the 
threshold for optical LO phonon emission. 

Optical and transport phenomena in QWRs have been of great interest lately. The great 
interest in transport phenomena JTB|-^TJ is due to the possibility that a 1DEG in a QWR will 
have the properties of a Luttinger liquid. The Fermi edge singularity |22| and the possibility 
of stimulated emission |23| enlarged interest in optics of QWRs . The important role played 
by LO phonons in the trapping of electrons |24| from three-dimensional extended states, 
into one-dimensional localized states, has been reported in Ref. |25| . 



Theoretical studies of QWRs were concentrated around elementary excitations of a 1DEG 
26| , |27[| , Fermi edge singularity 0], excitonic effects [p9H3~lfl, electron-phonon scattering rates 
[32| and interactions 0,0, relaxation of photoexcited carriers [j3~5| , calculation of envelope 
states P6fl , and conductance of QWRs with self-consistent broadening effects pTf . 



II. QUANTUM KINETIC EQUATIONS 



The retarded Green function and the "statistical" Green function which appear in the 
kinetic equations are defined as 



G r (r 1 ,a l7 t 1 ;r 2 ,a2,t 2 ) = -i9(t l -t2)( ^(ri,^,^),^^,^,^) ) 



(1) 



G s (ri, ci, t x \ r 2 , cr 2 , h) = -i{ *(r 1? en, *i), ^ f (r 2 , 2 , t 2 ) ) 



T 



(2) 



where \l/(r, 0, t) are the field operators in the Heisenberg representation. The square brackets 
with a plus (minus) sign signify an anticommutator (commutator). The upper signs refer to 
Fermions and the lower signs refer to Bosons (this convention is used throughout the paper) . 
The spin degrees of freedom a, will be suppressed from now on, and one can think of them 
as included in the space coordinates. 

In case of interacting particles the Green functions obey the Dyson equation 

G(r 1 ,£ 1 ;r 2 ,£ 2 )=G (0) (ri,£i;r 2 ,£ 2 ) + 

/ dr 3 dr 4 / dt 3 dt 4 G {0 \r 1 , t ± ; r 3 , t 3 )S(r 3 ,t 3 ; r 4 , £ 4 )G(r 4 , t 4 ; r 2 , t 2 ), (3) 

where G is a matrix of Green functions, and S is a matrix of self energy functions 0: 



G 



\ 



G a 
G r G s 



s a 



(4) 



/ 



The superscript (0) indicates free particle Green functions. The advanced Green function 
G a = G r *. 

In the case of free particles we can use the free particle eigenstates ip a (r), and the 
corresponding annihilation operators d a (t), to write \l/(r, t) = 'Y^ a %l) a {r)d a {t). In this case 
Gr^°)(ri, t\\ r 2 , t 2 ) depends only on the time difference t 2 — 1\, and can be Fourier transformed 
in time into G l ^°'(ri, r 2 ; e) (the index I stands for r, a or s). For an orthonormal set of states 
a we can define Ga '(e) such that 



G'(°>( ri , r 2 ; e) = J2 MriW a (r2)GT\e). 

a 

From the definitions (0), @ we obtain (near the singularities) 



(5) 



G r(0) { 



V{e-e a ) -ivr5(e-e Q ) 
sign(e) {V[e - sign(e)e Q ) - m5{t - sign(e)e a )} 



(6) 



and 



G s j°)( 



(7) 



—2ni(l — 2n a )5(e — e a ) 
—2m(l + 2n a )5(e — sign(e)e a " 

where e a is the energy of state a, and V is the principal part function. The upper (lower) 
term in the column corresponds to Fermions (Bosons). The occupation number of state a, 
n a = (o^da), is not necessarily the equilibrium occupation number. 

In the case of interacting particles in external fields, if one neglects coherence effects, the 
Green functions are still diagonal in a, but the energy levels are broadened and shifted. Let 
us first look at the time invariant case (systems under d.c. conditions). The level shift A, 
and the level broadening T, are defined as 



A a (e) = ReS^(e), T a (e) 



sign(ej 



ImE'fel 



(8) 



where J7 a is the retarded self energy of particles in state a. The column has the same 
meaning as in expression (y). 

The retarded Green function is an outcome of the Dyson equation (§), that is written as 

^( e ) = GL (0) ( e ) + GL (0) (e)S^(e)GL(e), (9) 

in the representation of the a states. Using the defintions @ we obtain from equation @ 
(for positive frequencies) 

G r a (e) = V(T a (e)\e -e a - A a (e)) - i7rS(T a (e)\e - e a - A a (e)), (10) 



where 



V{r a {e)\e-e a -A a {e)) = ta ^ ] ,, (11) 

1 Wl K>> (e-e Q -A Q ( e )) 2 + (r Q ( e )/2) 2 ' l ^ 



is the "smeared" and "shifted" principal value, and 

r Q (e)/27r 

( e -e Q --A Q (e)) 2 + (r„(r)/2) 

is the "smeared" and "shifted" delta function. 



S(T a (e)\e -e a - A a (e)) = A "\V 2 , /T1 , w ^ 2 , (12) 



Comparing expression ([H]) with expression (|6]) we see that G T a (e) is in fact a generaliza- 
tion of Gq (e), where the level shift and level width have been introduced into the principal 
value and delta functions. We therefore write G s a (e) as a generalization of Ga (e), thus 
defining the occupation functions n a (e) 

G B a (e) = -2m(l T 2n a (e))6(T a (e)\e - e a - A«(e)). (13) 

Expressions (|I^) and (^) were written for the case of positive frequencies, which will 
be assumed from now on. In case of negative Boson frequencies one has to multiply the 
retarded Green function and e a + A a (e) in these expressions by minus one. 

The definitions above agree with the relation between the statistical Green function and 
the retarded Green function at thermal equilibrium 

G s a (e) = 2z(l T 2n T (e))lmGl(e) } (14) 

where nr(e) is the equilibrium distribution function (Fermi or Bose-Einstein). The smeared 
delta function obeys the following normalization 

JdeS(r a (e)\e-e a -A a (e)) = l, (15) 

and the occupation numbers are obtained from the occupation functions in the following 



manner 



n a = / deS(r a (e)\e-e a - A a (e))n a (e). (16) 

Relations ( |I~4"D - ([UJ) can be obtained from the Lehmann representation of the Green 
functions. These relations are valid only if the one particle Green functions can be considered 
as diagonal in a, in other words, in the absence of coherence effects. 
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As a result of definitions @ and (jig), the three unknown functions ReG r , ImG r and iG s 
have been replaced by a different set of three unknown functions: r a (e), A a (e) and n a (e). 

If the system is not under d.c. conditions, one can define a new pair of time variables: 
the relative time £12 = £1 — £2, and the "center of mass" time £12 = (£1 + ^)/2. Only slow 
processes will be considered so that all quantities vary in £ with some time scale that is large 
compared to the inverse characteristic energy of the particles. In such a case one can still 
use the formalism given above for the Green functions, but with £ as a parameter of the 
problem, in addition to e (see for example Ref. 0). 

For slow processes 

Gl(t, e) = Gf\t, e) + Gf\i, e)S r Q (£, e)G£(f, e). (17) 

Again G T a (i, e) can be written as a generalization of Ga (e), in the same manner as before, 
but now £ is a parameter which appears in all functions. Thus, expressions (|T0|) and ([13]) 
for G T a and G s a , can be used with the three unknown functions T a (t, e), A a (£, e) and n a (i, e) 
that depend on £ as well as on e. 

The known procedure for obtaining a kinetic equation for n a (i, e) is to apply the operator 
G~2* — G~l to the Dyson equation for G s a (i, £12) 0, where the operator G~l is equal to 
id/dtt — e a for electrons, and d 2 /dtf + e 2 a for photons and phonons (we use the convention 
h = 1 throughout the paper). The slow variation in £ should be considered and then the 
equation should be Fourier transformed in the time difference coordinates. This procedure 
leads to the QKE 

^G s a (t, e) = -2£ s a (£, e) lmCF a (i, e) + 2ImS I a (£, e) G s a (t, e). (18) 

From now on the "bar" over £ will be dropped, with the understanding that all time variables 
are in fact center of mass time variables. Equation ( |T8"D for photons and phonons includes 
a factor of e/e a that multiplies the time derivative on the left hand side. Since we are 
interested in off-shell energies that are close to the on-shell energy, this factor can be taken 
to be one. 



In order to obtain equations for V a (t, e) and A a (t, e) one has to substitute an expression 
for the retarded self energy function, in terms of the one particle Green functions, into the 
definitions of T and A (|8|). These equations are in fact the imaginary and the real parts of 
equation (|iT|). This yields coupled, self-consistent equations for T and A, that are also n 
dependent. It is important to stress that the occupation function depends on time explicitly, 
while the dependence of V and A on time is only through their dependence on n. Therefore 
in order to find n one has to write a kinetic equation, that will include time evolution through 
time derivatives. The equations for the level shift and width will include their dependence 
on time only through the appearance of n, and will not include time derivatives. 

If one uses the self-consistent-Born approximation for the self energy functions, one finds 
that the equations obtained for n a (t, e), T a (t, e), and A a (t, e) could have been obtained easily 
by applying a few generalizations to the BE for n a . In the transfer from the BE to the QKE 
the number of unknowns increases from one (n), to three (n, T and A), therefore apart from 
the kinetic equation we require two more equations for the unknown functions T a {t,e) and 
A a (t, e), that appear in the equation for n a (t, e). 

We present here a recipe which allows one to go from the BE for n a (t), to the QKE for 
n a (t, e), and from the decay term in the Boltzmann equation (which is defined below) to an 
equation for T a (t, e). An equation for A a (t, e) can be obtained from the equation for T a (t, e) 
using Kramers-Kronig relations. 

The BE, in its most general form, can be written as 

d 

n a = ^2 l^a,/3,7.J 2 [(1 =F n a )np(l =F %)... - n a (l =F %))%...] 2n5(e a - ep + e 7 ...) = S. 



dt 

0,7- 



(19) 



This equation describes the evolution in time of the occupation number of particles in state 
a, due to interactions with particles in states (3, 7..., with matrix elements M a ^ n .... Different 
states can also mean different particles. The right hand side of the equation includes delta 
functions that are responsible for energy conservation. 

There are many cases in which it is convenient to think of the collision integral S as 
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composed of "scattering in" events and "scattering out" events: S = S{ n — S ont . The first 
term on the right hand side of equation ([19]) is the scattering in term: It describes processes 
in which particles enter state a due to the interaction with other particles. This term can 
be written as S{ n = (1 =F n a )Q a , where Q a includes all the terms on the right hand side of 
the equation that do not multiply n a . We refer to Q a as the "generation" term. The second 
term is the scattering out term: It describes processes in which particles leave state a. This 
term is usually written as S ont = n a /r a , where r a is the scattering out time from state a. 

In the context of the generalization of the BE to a QKE, it is more convenient to think 
of the collision integral S as made up of "generation" and "decay" terms: S = Q a — T a n a , 
where the decay term — T a n a includes all the terms containing n a . Thus r Q = ±Q a + l/r a 
is the total decay rate of particles in state a. It is the total decay rate, and not only the 
scattering out rate, that is related to the retarded self energy |Q. The general form for the 
decay rate T a , as it is deduced from equation fll9|) , is 



T a = J2 l M «,/3,7...! 2 i^M 1 =F <)••• + (1 =F np)n T ..\ 2vr5(e a - e p + e 7 ...). (20) 

0,7- 

In a d.c. situation one obtains for the occupation numbers 

n a = Q a /T a . (21) 

In order to transform the BE into a QKE one should apply the following rules. To the 
left hand side of equation fll9"D apply 

• n a(t) —> ( n a(t,t) T |) £(r a (£) e )l e — e Q — A a (t,e)), where e is the off-shell energy of 
particles in state a. 

To the right hand side of the equation apply the following steps 

• The occupation numbers for particles in states a, (3... are replaced with the occupation 
functions, that depend on the off-shell energies e, e'..., corresponding to the on-shell 
energies e a ,ep... respectively: n a (t) — > n a (t,e), ripit) — > ripit^e')... 
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• The summand under the sum over all states /3, 7... is multiplied by the product of the 
smeared delta functions for particles in these states with the corresponding off-shell 
energies e',e"..., and integrated over the off-shell energies: 

E "> E [de'de'\..6(rp(t,e%'-ee-A (t,e'))6(r^t,e")\e"-e^-A 1 (t,e")) x ... 

/3,7- 0,7- 

• The energy conservation delta function of the on-shell energies e a ,€p... is replaced with 
an energy conservation delta function of their corresponding off-shell energies e, e' ...: 
6(e a ...)-*6(e...). 

• Multiply the entire right hand side by S(T a (t, e)|e — e a — A a (t, e)). 

In order to obtain an equation for F a (t, e) one should apply the first three steps of the 
generalization rules for the right hand side of the equation, to the expression for r o (p0[). 
The QKE can be written generally in the form 

^ I (n a (t, e) t I) 8(T a (t, e)|e - e„ - A a (t, e))| = 

{£ Q (£, e) - T Q (t, e)n Q (t, e)} 5(r a (t, e)|e -e a - A a (t, e)). (22) 

The function Q a (t, e) contains all the terms on the right hand side of the QKE, that do 
not multiply n a (t,e). It is the generating term - the term in the kinetic equation that is 
responsible for the particle generation rate. Equation ( P2"D shows clearly that the particle 
level width T a (t,e) is in fact the total particle decay rate. Note that due to ( |i5|) the 
integration of equation ( ^2|) over e results in the disappearance of the =Fl/2 factor from the 
left hand side. The occupation function in the time independent case is given by 

n a (e)=g a (e)/T a (e). (23) 

The representation of the QKE in terms of n, T and A is more convenient than the 
representation in terms of Green functions, since one can clearly see all physical processes, 
and use the physical intuition that one gained from the BE in order to simplify the QKE. 
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FIG. 1. Hot luminescence from a QWR: Electrons are excited from the highest valence subband 
£fc, to the lowest conduction subband e/%, via photons v. These electrons relax to the bottom of the 
band emitting LO phonons lu . The hot luminescence v' is due to recombination of electrons from 
the vicinity of the threshold for phonon emission (dotted line) with holes in an impurity level — e 
(large dashed line). 
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III. LUMINESCENCE FROM A QUANTUM WIRE 

We now use the formalism given in section [TJ in order to deal with the specific problem 
of hot luminescence from a QWR. The case in which electrons are excited from the highest 
valence subband to the lowest conduction subband via photons, and then relax via LO 
phonon emission (see Fig. |I]) is considered. All other conduction and valence subbands are 
neglected (in contrast to Ref. Pj| , |2"%|j ). We consider a situation in which electrons are excited 
by a narrow band, noncoherent laser, just above the threshold for optical phonon emission 
€k = u , and describe the hot luminescence due to the recombination of electrons with holes 
in an impurity level. It is assumed there is some mechanism due to which the electrons leave 
the bottom of the conduction subband, so that there is no electron accumulation there, but 
this recombination can be neglected in the balance equations near the threshold. The spin 
degree of freedom will be completely disregarded. 

Unless specified otherwise we will assume that the electrons are excited above the thresh- 
old, since we are mainly interested in luminescence due to these electrons. Below the 
threshold the electrons relax to the bottom of the band by emitting acoustical phonons. 
The relaxation rate due to acoustical phonons r^ 1 , is much smaller than that due to optical 
phonons, and can be neglected above the threshold. 

The QWR runs along the z direction so that the electron wave functions are localized in 
the xy plane, and free waves in the z direction. The wave function in the conduction band 
can be written as e lkz (f) e (r)/L 1 ^' 2 , where L is the wire normalization length, and r is a vector 
in the xy plane. The electron energy is e^ = k 2 /2m e . All energies are measured from the 
bottom of the conduction band. The holes in the valence band have wave functions of the 
form e ipz cf) h {v)/L 1 ' 2 and energy e g + e p , with e p = p 2 /2rrih. 

The exciting photons are taken to be plane waves and the photon frequency is v — 
Vf = c|f| (f is the photon momentum and c is the light velocity). The LO phonons are 
three-dimensional and have a flat dispersion law, u; q = uo . 

The impurities are distributed randomly in the bulk. The wave function of a hole in 
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impurity i is of the form ip(r — Ti,z — Zj), where r«, Zi are the coordinates of the impurity 
position. 

In order to analyze the hot luminescence we neglect light polarization effects for simplic- 
ity, and assume that the wire and the crystal are cylindrically symmetric (an assumption 
that would not fit V-grooved wires nor the assumption of Ref. |56|] ). Thus, the wave functions 
of the emitted photons are 

X/, m ,n(r, z) = ( — 5 ] e^^U^nr), (24) 

where R is the normalizing radius and L the normalizing length of the crystal, and z, ip and r 
are the cylindrical polar coordinates. J m are the Bessel functions, and (K m ^ n R) is the n th zero 

1 /9 

of J m . The luminescence photon frequency is given by v' = v^ m ,n — c {f 2 + K m.«) (where 
the prime signifies this frequency is of luminescence, and not of excitation). We neglect 
all excitonic effects (these are treated for the case of optical absorption in one-dimensional 
semiconductors in Refs. |29|- |31|1 ). 

We begin by writing down the BEs for the electrons in the conduction band, and for 
the luminescence photons. It is shown why this description fails and one has to turn to the 
QKEs. We will then employ the generalization rules in order to obtain the QKEs from the 
BEs, and use the QKEs in the analysis of the luminescence. 

A. The Boltzmann Description 

We assume low temperature and a weak excitation level. At equilibrium, when there is no 
excitation, all the electrons are in the valence band, and there are no phonons (N(u ) = 0), 
due to the low temperature. In other words, there are no electrons in the conduction band 
(n(fc) =0), no holes in the valence band {n{p) = 0), and the impurity level is fully occupied 
by holes (n(e) = 1). 

The excitation creates electrons in the conduction band, and holes in the valence band. 
Most of the excited electrons relax to the bottom of the conduction band, emitting LO 



13 



phonons, but a small fraction of them recombine with holes in the impurity level, producing 
hot luminescence. By the assumption of weak excitation we mean that N(u ), n(k) and 
n(p) are small. The luminescence would be even weaker than the excitation, therefore n(e) 
remains close to one. Due to the weak excitation the luminescence is spontaneous, and the 
occupation numbers for the luminescence photons are also small. 

In the derivation of the BEs, terms that are nonlinear in the electron occupation numbers 
will be neglected. We will also neglect nonequilibrium contributions to the phonon and the 
hole (both in the valence band and in the impurity level) occupation numbers. 

The interaction of the electrons with the luminescence field is neglected in the balance 
equation for the electrons, in comparison with their interaction with the excitation field, due 
to the weakness of the excitation. In the case of c.w. excitation the balance equation for the 
electrons is 

= E E KT,f I' 2™>("f ~e k -e g - e p )(l - n(k))N(f) - 
f v 

E E l M SSI 2 2n6 ^ - ^ - <">(*)» ( 25 ) 

q k> 

where M exc is the matrix element for the electron-exciting-photon interaction, and M c ~ LO 
is the matrix element of the electron-LO phonon interaction. 

Let us first look at the electron-exciting-photon matrix element. In this interaction 
one photon of wave vector f is absorbed, and a hole of momentum p and an electron of 
momentum k are created. The interaction, in second quantized form, is given by 

I dv [dz^&Ur,z)A cxc (r,z)V h Ur,z). (26) 

J J m c 

The mass m is the bare electron mass and the constant p is the "bare" electron momen- 
tum operator (in the direction of the photon polarization, assumed linear for simplicity), 
sandwiched between the valence band top and the conduction band bottom Bloch wave 
functions. \l/ c (r, z) and ^/ h (r,z) are the electron and hole field operators, respectively. The 
photon field operator is given by 
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ir,z) 



£ 



2ttc 2 

T 7 ^ 



2\ 1/2 



e^+^Uf + alf 



(27) 



where fj_, / are the components of f in the xy plane, and in the z direction, respectively, df 
is the annihilation operator, and V is the normalization volume of the crystal in which the 
quantum wire is embedded. 

In order to find M cxc one should sandwich interaction ( p5|) between the states (n(k) = 
1, n{e p ) — 1, N{ — 0| and | n(k) = 0, n(e p ) — 0, Nt — 1) . This is equal to 

„2\ 1/2 



M^ f = ^- C^S-\ - [ dr /"dze- i(fc+p)z 0°*(r)0 h *(r)e^- r+l/z . 
m c \ Vvf J L J J 



(28) 



The photon momentum is small and so is r, since it is limited by the wire cross section, 
therefore e* f±r = 1. Integrating over the z coordinates the square of the matrix element is 



M exc 

\ Ivl k,p,f\ 



ep \ 27rc 2 



dr 6 e (r)6 h (r) 



S 



f,k+p- 



(29) 



m c) Vvf 

In the electron-LO phonon interaction one phonon of momentum q is emitted, an electron 
k is annihilated, and an electron k' is created. The interaction of the electrons with the 
polarization created by the LO phonons is given by 



f dr [ dz^^y e Ur,z)B(r,z)&(r, 

J J q 



*)> 



(30) 



where 7 is some interaction constant (see for example Ref. [3£|, |39|), and the phonon field 
operator is 

1 



1 ' ; ^ \2Vpuj ) 



1/2 



3 ?qx-r+J<j2 



a Q + a 



(31) 



The components of q are defined in the same manner as those of f , d q is the annihilation 
operator, and p is the reduced mass per unit cell in the lattice. 

The matrix element M e_LO can be found by sandwiching interaction (^) between the 

states (n(k) = 0,n(k') = 1,N(lu ) = 1| and \n(k) = l,n(k') = 0,N(u o ) = 0). Carrying out 
the integration over the z coordinates, the square of the matrix element is 



| A/f e-LO|2_ 1 4vraa; 

I fc,fe',q I 



3/2 



\q\ 2 V{2m c y/ 2 



drl^irfe-^- 1 



4 



k—k',qi 



(32) 
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where we exchanged 7 for the known Frolich constant a, using the relation 7 2 = 
(p Wo 2 /47r)(2^ /m c ) 1 /2( a / e 2) (see Ref. [0). 

Substituting |M exc | and |M e_LO | into equation (p5|), we can sum over p and q. Taking 
the normalizing volumes to infinity and thus exchanging the sums with integrals we obtain 

° = IwT \M c * c (vf)\ 2 2n5(v f - e k -s g - e/ _ fc )(l - n(k))N(f) - 

a\M c - LO (k-k')\ 2 27r5(e k -e k ,-cu )n(k), (33) 

J In 



where 



|M» W | 2 = ^ 

m;z/ f 



dr<f(r)0 h (r) 



(34) 



f c-LO- -mi2 -1 47ro;a; 



3/2 



|M—(* -*-)! =_p_^^-n, (35) 



and 



a(k-k')= J dt f dr'\4>°(r)\ 2 \(f> e {r')\ 2 K (\k-k'\\r-r'\). (36) 

The zeroth order Bessel function Kq is a result of the integration over q^. 

Since the matrix elements fl34|) and ( p5|) are smooth functions of their arguments, and 
we are interested only in processes which involve electrons that were excited close to the 
threshold, we may take them at the threshold values. The threshold value of v? is v — 
e g + (1 + f])uJo, where rj = m^jm^. This is the frequency that will excite electrons from the 
valence band to the threshold in the conduction band. The threshold value of \k — k'\ is 
q = {2m c uj ) 1 / 2 . 

The electron wave functions that appear in a limit the spatial integration to r < a, 
where a is the wire width, therefore we can estimate that <j(q ) oc | ln(g o)| for q a <^C 1 and 
<j(q ) oc l/(q a) 2 for q a ^> 1. For GaAs q = 2.5 x 10 6 cm _1 , thus for most wire widths we 
are in the regime of q a 3> 1, in which a is expected to decrease with increasing wire cross 
section. The constant aa(q ) = a* is the effective Frohlich constant in the one-dimensional 
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case. According to Ref. |32|] a* ~ 0.1 (which is of the same order as the bulk Frohlich 



constant) for a wire of cross section 200 x 100A , while according to Ref. [^] a* « 0.02 for 
a wire of cross section 300 x 100 A . 

We now perform the integrations which appear in equation (^). We begin with the 
generating term Q(e k ), the free term (one that doesn't include n(k)) in the balance equation 
(J33J) for the electron occupation number. This term is physically the generation rate of 
electrons in the conduction band. Since / <C k we can approximate £/_& ~ £_&. Using 
df = v\dvid£lij(? « v 2 dvt dfif/c 3 , where dfif stands for the solid angle increment, we 
obtain 

^ \M^{v)\ 2 2n5(v { -e k -e g - e- h )N(f) = 

2tt\M cxc (u)\ 2 D(u)(N(e k + e g + 6^ k )) =G(e k ), (37) 

where D{y) = z/ 2 /(27r 2 c 3 ), is the photon density of states, and 



(N(e k + e g + £_*)) = J ^iV(f) 



(38) 

v f— e fe+£g+£-fe 



is the angular average of iV(f), for photons of frequency e k + e g + E- k , that create electrons 
with energy e k . 

We consider a narrow band photon excitation field such that the angular average of the 
occupation number is given by 

wo>=j(g)' 7 y, ,,, , (39) 

where J is the excitation field energy density, Az/ is the spectral width, and v is the central 
frequency. We are interested in a narrow band excitation close to the threshold, therefore 
Av <^C u) , and the detuning (the difference between the central excitation frequency and 
the threshold frequency) is small: v = v — v <^ u . 

The integrations that appear in the expression for the decay rate T(e k ), the term that 
includes all coefficients of —n(k) in the kinetic equation ([33]), will now be carried out. The 
first contribution to T(e k ) is that of the photons, and is equal to the generating term (P7|). 
The second contribution to r(ejfc) is due to LO phonons, 
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1 /9 

r L o(e fe )= [^\M e - LO (q )\ 2 2n6(e k -e k ,-u; ) = 2a*u; (^^-) . (40) 

J *7T \e k — IX> J 

The photon contribution to the decay rate is negligible compared to the phonon con- 
tribution since radiative processes are slow compared to non-radiative processes. Therefore 
r(efc) can be written as 

r 3 \ 1/2 



r(e fc ) = r LO (e fc ) = c — , (41) 

where 

T c = {2a*) 2 ^ . (42) 

The meaning of the energy scale T c will be explained at the end of this section. For a* = 0.02 



1 35| one obtains T c = 4meV, thus T c ^ lo = 36meV. 

Below the threshold, where IYo = 0, the acoustical phonon contribution to T is impor- 
tant. Since the latter is a smooth function near the threshold, a good approximation for T 
is T(efc) = r^ 1 = const at e^ < u . 

The occupation number of the electrons that were excited above the threshold is given 
by expression (pl|), where Q and T are given by expressions ( pTj) and (|4T| ) respectively. It is 
clear from the expressions above that n(k) depends on k through e& only. 

We now turn our attention to the kinetic equation for the luminescence photons, 

-N{f,m,n) = J2 W l ^ m , n \ 2 2ix8{v f , m , n - e k - e) x 

k 

[(1 + N(f, m, n))n{e k )n{e) - N(f, m, n)(l - n{e k )){l - n{e))} , (43) 

which for n(e) = 1 and in case of spontaneous luminescence is 

-N(f,m,n) = J2 \M{j m j27r6(u f , m , n - e k -e)n{e k ). (44) 

k 

M lum is the matrix element for the luminescence photon-electron interaction. The time 

derivative in the equation above is kept in order to clarify what is the luminescence source. 

Let us first look at the matrix element that describes the recombination of an electron 

with a hole from a specific impurity. We will then have to sum over all impurities and average 



over all impurity configurations. In this process a photon of quantum numbers /, m, n is 
emitted, and an electron k and a hole in impurity i are destroyed. The interaction is given 
by 



d r I c^^^ e (r,z)i lum (r,;z)^ h (r,z) 
m c 



where 



/ 2vrc 2 \ l l 2 / \ 

i lum ( r , z) = 22 [ (x/,m,n(P, Z)af,m,n + X/,m,n( r , z Y a \,m,n) 

j,m,n v J ' ' ' 

is the field operator of the luminescence photons. 

The matrix element of the interaction above between the states (N(f, m, n) 
0,n(e) = 0| and \N(f,m,n) = 0, n(efc) = l,n(e) = 1) for holes in impurity i is 

27TC 2 



(45) 



(46) 



l,n(e k ) 



1/2 



MkJ,m,nV> 



ep 



m c y<KR 2 L 2 v fjmin \J m+1 (K min R)\' 
dr I dze tkz ct> e {Y)i,{v-Y %1 z- Zl )e- tfz ~ imip J^ mtn r). 



(47) 



We now sum the matrix element squared over all impurities and average over all impurity 



configurations. This yields 



M 



lum 
k,f,m,n\ 



ep 



27TC 2 



X 



m cj TiR 2 L 2 v f ^ n \J m+1 {K m ^R)\" 

^vl dTl [ dZl [ dT [ dz eifc2 ^( r M r -n,z- Zi)e- lfz - imv r m (K mtn r) 



(48) 



Carrying out the integrations over the z coordinates we obtain 

2 r. 2 



M 



lum I ^ 



k,f,m,n I 



ep 



m r 



2ttc z 



W< 



o c J irR L lSf } m,n\Jm+l{ K m,nR)\ 



dr c (r)^(r - r 4 , -k)e~ im * r m {K m>n r) 



T L x 



(49) 



where ip{r,k) = J dze ip(r, z). The second argument of ip, f — k was replaced by —k 
because f <^ k. 

Since the impurities are randomly distributed in the bulk, the average over all impurity 
configurations of the product of the hole wave functions, that appear in fl4"9"|), can depend 
only on the difference between their coordinates. Thus we define 
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£(r- r',-k) = / driip(r - ri,-k)ip*(r' - Ti,k). 

The summation over impurities divided by the normalizing volume renders a factor of 7ii mp 
- the impurity spatial density. The averaged matrix element squared can then be written as 

| M lum |2 = (ePo\ 2 2^ f j 

1 fc ' / ' m '" 1 V"W nR*Lu fjm , n \J m+1 (K m , n R)\ 2 imp J J 

c (r)0 c *(r')/:(r - r', -k)e~ im ^-^ J^(K m , n r) J m (K m , n r'). (50) 

The matrix element squared is much larger for m = than for other values of m, since r 
and r' are constrained to a small region in the xy plane (due to the electron wave functions), 
and Jo is the only Bessel function which is finite at r — > 0, Jo( r )| r _>o = 1- Thus the matrix 
element squared is proportional to 5 m ,o- From now on the index m is omitted with the 
understanding that we are dealing only with m = 0. 

The matrix element squared is a smooth function of v^ n and k, therefore these can be 
substituted by their threshold values. The threshold values of i>f >n and k are v' = e + u and 
k = {2m c uj ) 1 / 2 respectively. The kinetic equation for the luminescence photons can then 
be written as 

^N(f,n) = J —\M l ™(k )\ 2 2n5(v f , n -e k -e)n(e k ), (51) 

where the matrix element squared is 

|M^ m (A; )| 2 = f^-X — ,n imp Lx 

1 n V ° n \m cj KB?W\J x (K n R)\ 2 P 

( dt I ' dY'<jf{r)<f?*{Y')C{r-r',-k ). (52) 

Substituting n(e^) into equation ( |5"I| ) and integrating over k we obtain 

|iV(/, „) = \M'r(K)f (-^-) ^ p^\- (53) 

dt \ u f,n-eJ r(i// in -£) 

The luminescence source can be characterized by the spectral dependence of E(v')dv', 
the energy of the emitted field, of all spectral modes within the interval dv' ', per unit length 
of the wire, per unit time. Multiplying the generating term of the luminescence photons 
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(the right hand side of equation (|53|)) by 5(u' — Vj^v'dv' / L and summing over / and n, we 
obtain 

ew) rf,' = ^EE K"(«r (1^7) " 2 ^7^ *' - "/-) **■ < 54 > 

The variable v^ n in the square root can be replaced with its threshold value e + u , due 
to the smoothness of this function. Taking the normalization volume (that appears in the 
matrix element) to infinity the sums over / and n are transformed into integrals. Performing 
the integration we obtain the final result 

E{y')=Cv'D{y')i{v'), (55) 

where 

c =fe) ^n bBp JdrJ<h'<ir{r)<r{T r )£{r-T',-k ), (56) 



is a constant, and 



' C.fn 1 r-\ / Orrr, \ V^ Q— - I /I /f6XC / f" , \ I 2 



2m e V' Q{y' - e) (2m c y /z 2tt \M exc (u)\ z D{p) 



e(i/) = | ^ j ^__j_ = | ^e j j ^_pL_^ Z(iV(( i H . 77) ^_ £)H _ £g))) 



(57) 



is essentially the product of the electron generation rate and the electron life-time at the 
corresponding energy. 

Since u'D(i/') is a smooth function of u', the spectral distribution of the luminescence 
is given by £(z/). Let us first consider a narrow band excitation that does not overlap the 
threshold, i.e. v ^> Av . In this case £(f') is a peaked function of z/, reproducing the shape 
of the excitation. It is of a Lorentzian shape of width Az/ = Az/ /(l + rj), centered at v' + u' Q , 
where u' = v /{\ + rf) (see Fig. ||) 



Vo V /2 Ai/ /(2tt) 

:i + ^) 3 r c V (z/ - ix - u' f + (Az/ /2) s 



£(" ) <* ( 7i I „^™ ) ~ ~ - M 2 , , A 77^2 - ( 58 ) 



The rescaling of z/ compared to z/ follows from the obvious relations (see Fig. [I]) e& + £& 
z/ — e g and z/ = e + e&. 
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v — V 



u'-u' 



FIG. 2. The behavior of £,{v') as a function of v' for v 3> Av (rj = 1/3 and v = 8Af ), where 
v' a = z? /(l + 77) and A.v' = Af /(1 + rj). An enlargement of its behavior for v ' — v' -C z>„ is shown 
in the inset. 




Ai/ /2 

FIG. 3. The behavior of £,(v') as a function of v' for Az^ D 3> t'o (rj = 1/3 and Af D = 8v ), where 

Au' = Au /(l+ V ). 

Only close to the threshold, when v' — v' <C £ /(l + ??)> the spectral behavior of the 
luminescence differs from that of the excitation and is given by (see the inset of Fig. 0) 






(59) 



When the narrow band excitation overlaps the threshold, i.e. v <C Au , the spectral 
distribution of the luminescence is very different from that of the excitation (see Fig. ^|): 

2 



£(*/) ex 



nAu T 3 J 2 



[y'-V') 1 ' 2 when */-P'<^°. 



(60) 



and 
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£(z/)oc AU ° ^A v-v')-^ 2 when v' - v' > ^. (61) 

2 7 r(l + r / )2rf 2V 2 

If one is interested in luminescence due to electrons that were excited close to but below 
the threshold, via an excitation that is centered above the threshold, i.e. v' — v' < and 
u > 0, then T(i/ — e) in expression (|57| ) should be replaced by r" 1 . In this case £(^') 

\3f3^V2 



has the same Lorentzian form as (|58|) with the amplitude (z/ /((l + r/) 1^)) replaced by 

r ac /(l + v)- 

One can expect that the luminescence distribution £(V)> derived from the BE, is not 
correct for frequencies v' close to the threshold v' . Luminescence of such frequencies is due 
to electrons with energies e^ near the threshold, that have a large width T(efc), (see equation 
(PI)), while the BE assumes that the width of the states is small compared to other energy 
scales of the problem. In our case one of these energy scales is e^ — uj - Equating the width 
with this energy scale, e^ — uj = T(efc), one finds an additional energy scale of the problem 
(and a critical decay rate), that is given by T c that has been defined before ((4^). 

When ek — uj < T c , we have T(tk) > e^ — u , and the assumption of the BE breaks down. 
One would therefore expect that for e& — u < T c in the description of the excited electrons, 
and for v' — v' < T c in the description of the hot luminescence, the results predicted by the 
Boltzmann equation would fail. We will learn from the QKEs that in fact the situation is 
more complicated. 

B. The Quantum Description 

As we have seen in the previous section, the Boltzmann description may fail for electrons 
excited to the vicinity of the threshold and hence for luminescence photons emitted by these 
electrons. Therefore, close to the threshold, the electrons and the luminescence photons 
have to be treated using QKEs. These QKEs will be written by applying the generalization 
rules to equations (|33D and (|5lD . It is assumed that the energies of the holes, the phonons 
and the excitation photons are not broadened. 
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In the generalization to the QKEs only those terms that were kept in the Boltzmann 
description will be retained. Threshold values of the arguments of the matrix elements will 
be substituted as was done in the Boltzmann description, due to the smoothness of the 
matrix elements. 

The QKE for the electrons is 

f^jde'5(T(k',e')\e' - e k , - Ae(k',e')) \M C - L0 (q )\ 2 2n5(e - e' - u )n(k,e), (62) 

where e is the electron off-shell energy, n is the electron occupation function, and F and Ae 
are the electron energy level width and shift, respectively. 

The coefficients of —n(k, e) in equation (|62"D can be recognized as the electron level width 
due to LO phonons 

r LO (M = J^ J de'5{Y{k'^)\e' - e k , - Ae(k',e')) |M e - L0 (g )| 2 2n5(e - e' - u ), (63) 

and the total electron width is given by 

T(k, e) = r LO + 1/Tac- (64) 

The photon contribution to the level width was neglected (as it was neglected in the Boltz- 
mann description of the decay term). 

One can see that the generating term in equation ( |62|) is equal to expression fl37|), the 
generating term of the electrons' BE, with the on-shell energy e k replaced by the off-shell 
energy e (in the exciting photon occupation number) 



g(k,e)=2n\M cxc (u)\ z D(u)(N(e + e g +e. k )). (65) 

The occupation function n(k, e) is found by substituting Q and Y into (^3|). The generat- 
ing term Q(k, e) depends on k through e_ k = r\t k only, therefore Q(k, e) = G(e k , e). It follows 
from ( |63D that T is independent of k (and so is the level shift). Therefore, n(k, e) = n(ek, e). 
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We) 




e- uj 

— r — r r r 

FIG. 4. The behavior of the electron level width due to optical phonons, at the vicinity of the 
threshold for LO phonon emission, when the electron level width at the bottom of the subband is 
considered. 



Equation fl63|) is a self consistent equation for the level width. Integrating the right hand 
side of the equation we find that 



We) 



a uj 



3/2 



e — u — Ae(e - u ) + ~r(e - u ) 



-1/2 



+ 



UJn 



Ae(e — uj n 



-Vie — uj 
2 v 



-1/2' 



(66) 



When e-w Cw„ the arguments of T(e — uj ) and Ae(e — u ) are close to the bottom of the 
band. In this case the level shift can be neglected since it is just a small renormalization. The 
level width close to the bottom of the band T , is due to scattering with acoustic phonons 
and due to thermal recombination. The latter contribution to T can be neglected, while 
the acoustic phonon contribution can be considered as a constant (though not the same as 
the width due to acoustic phonons at the threshold). We assume that T is small compared 
to all other energy scales of the problem. Thus close to the threshold 

r LO (e) = «x 3/2 



e - UJn 



^ 



-1/2 



+ 



UJn 



^ 



-1/2' 



(67) 



The behavior of T-^o(e) when e is close to uj is shown in Fig. f|. 

From expression ( |67|) we see that above the threshold for e — u> ~^> r o /2, T(e) can be 

written as 

1/2 / r 3 \ 1/2 



r(e) = r LO (e) = 2a*uj 



UJn 



e-uj 



UJo 



(68) 
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Comparing expression (|68|) to expression (|41|), it is evident that above the threshold T(e) is 
equal to T(efc) when the on-shell energy e& is replaced by the off-shell energy e. Since T is 
the smallest energy scale in the problem, expression (^) can be used for the electron level 
width above the threshold. 

Contrary to the Boltzmann description, the electron-LO phonon scattering contributes 
to the electron level width below the threshold as well. When uo — e 3> T the electron 
width due to optical phonons decays like |e — u \~ 3 ^ 2 , much faster than the decay above 
the threshold (see Fig. f|). The contribution of the acoustic phonons to the electron width 
below the threshold r^ 1 , will be negligible compared to IYo for uj — e <C (r o r ac ) 2//3 r c /2. 
Since r o r ac <C 1 and r c r ac ^> 1 the right hand side of the inequality above is much smaller 
than T c , and much larger than T (see the crossing point of r L o(e) and r^ 1 in Fig. [|). 

The QKE that is obtained by applying the generalization rules to the BE for the lumi- 
nescence photons (|5lD is 
d f/1 



()f [K2 + N(f, n, v')\ 8{<y(f, n, u')\u' - v f , n - Au(f, n, u')) 

KlUi n > v ')\ v ' - v f,n - Ai/(/, n, v')) X 

J^Jde \M l r{k ) | 2 S(T(e) \e - e k - Ae(e))2n6(u' - e - e)n(e k , e), (69) 

where v' is the photon off-shell energy, N is the photon occupation function, and 7 and 
Au are the photon energy level width and shift, respectively. All the photon functions are 
written for m — 0. 

By applying the generalization rules to the right hand side of the full BE ( f43| ) one can 
find the photon level width 

l(n,v') = -J^Jde\M i r(k )\ 2 5(T(e)\e-e k -Ae(e))x 

2ic6(i/-e-e)n(e k ,e). (70) 

Note that 7 < corresponds to photon generation. It is evident from expression (|T0"| ) that 
7 does not depend on /, and therefore neither does Au. 
Using (|70|) equation ( |69|) can be written as 

26 



-7(71, v')6("/(n, v')\v' - Uf tTl - Av(n, v')). (71) 

The right hand side of equation ([71]) is the photon generating term, the term from which 
the luminescence spectral distribution will be obtained. 

As in the classical case, the luminescence source can be characterized by the spectral 
dependence of E(v')dv' ', that is obtained by multiplying the right hand side of equation ( |7T| ) 
by v'dv' jL and summing over / and n 

E(v') dv' = -j J2J2 7 ( n ' v'M'yfa v ')\ v> ~ v fv ~ Az/ ( n > z/ ')) dv '- ( 72 ) 

n f 

The main contribution to E(y') comes from / and n such that v^ n is close to v' . 
It is convenient to write j(n, u') in the following form 

7(ny) = - m , T \ m a C ^% (73) 

where C is given by (^6[), and 

ZW = j^j de 6 ^ e )\ e - e * - A e ( e ))27r<5(^ - e - s)n(e k , e). (74) 

When the electrons' energy is taken to be on shell, that is the electron smeared delta function 



is reduced to a singular delta function, £(z/') is reduced to expression (|57|) . In writing 7 in 
the form (|73|), the normalization volume has been separated from the rest of the function, 
that is independent of the index n, and can therefore be taken out of the summation in 
E(v>). 

Taking the normalization volume to infinity the sums over / and n are transformed into 
integrals. From equation ( |73D it is evident that 7!^^^ —> 0. This is due to the fact that 
the photons are emitted into an "infinite" space, while their interaction with the electrons 
is confined to the finite volume of the wire. In this case 7 and Au are negligible compared 
to the other widths, and the smeared photon delta function that appears in the expression 
for E(v') becomes a singular delta function, giving 
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/°° fif f°° k cIk 
^l ^b{y<-v Ln ) = Cv<D{yi)i{y'). (75) 

Since v'D(v') is a smooth function it is clear at this point that in order to understand the 
dependence of E(i/') on z/, the behavior of £(V) should be analyzed. 

Performing the integration over e in (0), and substituting the expression we got for 
n(e k ,e), we get 

, A , , U9 2n\M cxc (u)\ 2 D(u) f°° de k .,„, Nl , 

(JV(i/-e + e fc + e g )>. (76) 

Since e k — uj <C u , we substitute oj for e k in the one-dimensional density of states. This 
integral contains in fact two Lorentzians. The first is of width T(i/' — e): The electron level 
width at energy e = v' — e. The second Lorentzian is of width Av /r]\ A width that is 
proportional to the spectral width of the exciting photon source. Since the excitation is 
such that e k is close to u> , and all the energy scales that can characterize the level width 
are much smaller than u , the integrand goes to zero when e k approaches zero from above. 
Thus, the lower boundary of the integral can be taken to minus infinity. As a result the 
spectral distribution of the luminescence is given by the function 
_ (2m c \ l/2 2vr|M-Q7)| 2 D{p) / C x3 

de ( ry-£)/27r \ / Az^vr \ 

, k \{v'-e-e k Y + T{v'-ef/A) \[(y> - e + r}6 k + e s ) -v \* + M*/±) ' l J 

The region of validity of the BE is evident from equation flT7|). In order to obtain 
the result of the BE ([57]), one has to replace the first Lorentzian with a delta function. 
This can be done only when its width is smaller than that of the second Lorentzian, i.e. 
r(z/ — e) <C Au /rj. It follows from this inequality that the Boltzmann description of the 
luminescence spectra is correct only far from the threshold, when v' — v' ^> r] 2 Y\j Av 2 = v c . 
The quantum interval u c differs from the naive estimate T c (see the discussion at the end 
of subsection FCI A[): It is larger for "narrow" band excitation (Az/ <C T c ), and smaller for 



"wide" band excitation (Az/ ^> T c ). 
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For an excitation that is mostly of frequencies above v + u c , the spectral dependence 
of the energy of the emitted photons will behave according to the predictions of the BE 
(see Figs. @-|3|). However, as long as the excitation is above the threshold, very close to the 
threshold £(V) will increase linearly with z/, and not as a square root, see (|59|). 

Non-Boltzmann behavior is obtained when Y(v' — e) ^> Au /i], that is v' — v' <C z/ c , and 
the second Lorentzian can be treated as a delta function. In this case one obtains 

n -1 

(78) 



CO') oc 



2n(l + r]) 2 



■/-P--A-V+ " 2 r? 



1 + 77/ 4(1 + r/) 2 v 1 -v> 

We bring here two specific examples of extreme non-Boltzmann behavior of ^{y'\ In both 
examples the excitation is centered within the quantum interval, i.e. v <C u c , and "narrow" 
band, i.e. Au <C T c . Due to the latter inequality v c ^> T c . 

In the first case the detuning is large, v ^> T c , therefore the energy width of an electron 
excited by the central excitation frequency, e^ = oj + v j{\ + r/), that is of the order of 
(rj?/z> ) 1//2 , is small compared to its distance from the threshold. This is true for most of 
the electrons. In such a case the electron states are "well defined", but the prediction of the 
BE for the electron distribution is wrong, since its width Az/ /(l + r\) is much smaller than 
the width of the states. As a result it follows from (^) that the luminescence spectra is 
symmetric and centered at the classical position v' + i>' , with v' Q = v /(l + rj), but its width 
is given by the quantum-mechanical width of the electron states (see Fig. |5|) 



Au' = ^L 



n 



W0- + V) 



1/2 



In the second case the detuning is small, v -C r c , and the electrons are excited to 
"badly defined" states, since the broadening of these states is larger than their distance 
from the threshold. As a result the luminescence spectral distribution differs greatly from 
that predicted by the BE. One finds from ([T8|) that (see Fig. ^]) 

2 
£(z/) oc — — - iy' — v) when v' — v' <S r c , 
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FIG. 5. The behavior of £(V) as a function of v' for v c S> v Q 3> T c 3> Ai/ (7/ = 1/3, P = 15Af , 
and T c = 5A^ D ), where v' Q = v j(\ + 77) and Az^ = 77 iT^,/v' ) 1 ' 2 /{I + 7/). An enlargement of its 
behavior for v' — v' <C (t/ 2 / 4)1^ / 'f/^ is shown in the inset. 




v — v 



FIG. 6. The behavior of £(^') as a function of 1/ for v c 3> T c ^> z> , Az/ D (7/ = 1/3, P = 2Az/ , 
and T c = 15Az> ). 
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and 

^ ,)cc 2n(i\r ] y vhy when »c»v'-»'»rc- 

In conclusion let us note that the hole dispersion plays an essential role in determining 
the quantum behavior of the luminescence spectra. The width of the second Lorentzian in 
(\n\) Auo/rj, is not equal to the width of the classical luminescence line Az/ /(l +r)), and the 
integral is not a simple convolution of the classical luminescence profile with the spectral 



function of an electron state. One can see from ([77]) that for a flat hole band (77 = 0) the 
second Lorenzian does not depend on the integration variable and the first Lorenzian is 
integrated to one, restoring the BE result, unexpectedly. 

The luminescence due to electrons that were excited below the threshold (e < uo ) by an 
excitation centered above the threshold (z> > 0) is not obtained by the simple exchange of 
Tlo by r^ 1 , as in the classical case. As long as v' — v' < (r o r ac ) 2 / 3 r c /2 one should take 
T = r LO from expression ( |67D since T LO (v' — s) > r^ 1 , and therefore although the electrons 
are below the threshold, it is the optical phonon contribution to the electron width that is 
dominant. 

IV. CONCLUSIONS 

We have presented a recipe that allows one to generalize the BE to a QKE. The QKE 
obtained by employing this recipe is the same as that obtained by the Keldysh Green function 
technique in the self-consistent Born approximation. The advantage of this method for 
writing the QKE is that it provides a physical understanding of the terms in the QKE, 
and it allows one to neglect those terms that were negligible in the Boltzmann description. 
This is due to the fact that the equations are written for quantities that are similar to the 
quantities described by the Boltzmann equation. 

We considered the specific example of hot luminescence from a QWR. The recipe de- 
scribed above was used in order to generalize the set of BEs that describe the problem to 
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a set of QKEs. Solving these equations we were able to describe the luminescence spectral 
distribution. 

We have shown that there is a domain of luminescence frequencies, that correspond to a 
domain of phtoexcited electron energies, for which the quantum description of the lumines- 
cence spectral distribution leads to a different behavior than that given by the Boltzmann 
description. This quantum domain could not be easily guessed from level width considera- 
tions. Two other nontrivial conclusions were obtained. The first is the role played by the 
hole mass in the definition of the quantum domain. When the hole dispersion relation is 
flat the quantum domain shrinks to zero and the results of the Boltzmann description are 
retrieved. The second nontrivial conclusion is that there is an energy domain below the 
threshold for LO phonon emission, in which the LO phonon contribution to the electron 
level width is dominant. 
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